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Abstract 



We define, in a consistent way, non-local pseudo-differential opera- 
tors acting on a space of analytic functionals. These operators include 
the fractional derivative case. In this context we show how to solve 
homogeneous and inhomogeneous equations associated with these op- 
erators. We also extend the formalism to d-dimensional space-time 
solving, in particular, the fractional Wave and Klein-Gordon equa- 
tions. 



2 



1 Introduction 



Interest in non-local field theories has always been present in theoretical 
physics, associated to several different motivations. 

In Ref. J. A. Wheeler and R. P. Feynman considered a description 
of the interaction between charged particles, where the electromagnetic field 
does not appear as a dynamical variable (action at distance). 

More recently, efforts were made to use non-local theories in connection 
with the understanding of quark confinement and anomalies and in 

string theories containing non-local vertices Ql^. 

Besides the possibility of non-local interactions, a field theory can also 
display non-local kinetic terms. 

Before renormalization theory became well established the possibility was 
considered of formulating finite theories by means of non-local kinetic la- 
grangians. A. Pais and G. E. Uhlembeck ^ were one of the first to analyze 
non-local theories in this context. 

The analytic regularization method introduced in Ref. 0, can be thought 
of as associated to non-local kinetic terms in the lagrangian, specifically to 
fractional wave and Klein-Gordon equations. 

This type of non-locality also arises in effective theories, when integrating 



over some degrees of freedom in an underlying local field theory |T^ . 

At a classical level, non-local equations containing arbitrary powers of 
the D'alambertian (□) have been studied in Ref. In this reference a 



non-trivial relation between the number of dimensions and the power of the 
operator was established to satisfy the Huygens principle (see Ref. [0). In 
particular, in 2 + 1 dimensions the usual wave equation leads to a Green func- 
tion that does not satisfy Huygens principle, while a non-local equation with 
□ does satisfy this principle (see also Ref. |13[)- It is not by chance that 



the pseudo-differential operator D^/^ appears also in the process of bozoniza- 
tion in 2 + 1 dimensions. In Refs. and |15] a mapping was established 
between Dirac's kinetic term and bosonic terms involving D^/^. In Ref. |]T^ it 
is shown how similar terms appear when 3 + 1 dimensional QED is projected 
to a physical plane. The kinetic term contains D"^/^, which in the static 
limit reproduces correctly the r^^ Coulomb potential, instead of the usual 



logarithmic behavior of 2 + 1 QED. This fact was first noticed in Ref. |12]. 

Also, in Ref. [|17|, a fractional generalized Fokker-Planck-Kolmogorov 
equation is proposed to describe anomalous transport in Hamiltonian sys- 
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terns, and in Refs. [T^l^ some particular Green functions for fractional 



diffusion and fractional wave equations were obtained. 

In Refs. []T3[ and , solutions to some particular non-local homogeneous 
equations were proposed on physical grounds. 

In this paper, we present a general mathematical approach to deal with 
pseudo-differential operators which enables us to obtain the whole space of 
solutions to homogeneous equations associated to them. Therefore, for a 
given pseudo-differential operator, it is possible to define a general Green 
function (having a general prescription to avoid the singularities). 

In section 0, we define the functional spaces in which we build-up the 
framework for a correct description of pseudo-differential operators. In these 
spaces we introduce the representation of analytic functionals as "ultra dis- 
tributions" p2| , a convenient way to handle and operate with the usual 
Green functions of Quantum Field Theories. 

In Section ^ we define local and non-local operators acting on ultra dis- 
tributions. We write and solve formally, homogeneous and inhomogeneous 
equations associated to them. 

In section H, some examples are examined, such as rational and power 
functions of the derivative operator. Solutions of the respective homogeneous 
equations are given. Also, Green functions are found, and their relations to 
fractional derivatives of Ref. are exhibited. 

In section ^, we extend our developments to space-time, by defining non- 
local functions of □. Also, fractional wave and Klein-Gordon equations are 
introduced, and solutions for the respective homogeneous equations are given. 
Corresponding Green functions are found. 

Finally, section ^ is devoted to a discussion of the developments and 
results of the paper. 



2 Functional Spaces 

We will start with the space ( of entire analytic test functions (p{k), rapidly 
decreasing in any horizontal band. We will call "ultra-analytic" , any function 
<^{k) G C- They are Fourier transforms of the space ( of all C°° functions 
(p{x), such that exp{q\x\DPip{x)) is bounded in R for any q and p. 

In view of the latter property, ( D Z, where Z is the space of Fourier 
transforms of K {C°° functions on a compact set). ([^, ch. 2, §1.1) 
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The dual of ( is the space (' of hnear functionals defined on (. In (' we can 
represent the propagators of a quantum field theory as analytic functionals, 
with the physical properties that are expected from them 

The general form of an analytic functional is: 

^((^) = f tPikMk) dk (1) 

The "density" ipik) is an analytic function and the line L can be deformed 
as long as no singularity of ipik) is crossed. L can be an open line or a closed 
loop. 

Not only L can be deformed without altering ipif)- Also, the structure 
of singularities of the density can be modified. For example, when ip{k) = 
(A; — r)~\ the development: (fc — r)^^ = J2n t"A;~"~^, allows the pole at = r 
to be represented by a series of multipoles located at the origin. This fact is 
closely related to the expansion of the analytic functional 6{k — r) as a series 
of (H, ch. 2, §1.4 ). We can see then, that an analytic functional 



can be expressed in more than one way. 



We are going to use, systematically, the following representation ||2T[l [p^] : 



^Piip) = dk mm e C' , ^ G (2) 

where ip{k) is analytic in {/c G C : \Im k\ > p} and 'ip{k)/kP is bounded 
continuous in {/c € C : \Im k\ > p}, p depending on ip{k), p & N {N = set 
of entire numbers). 

The path F runs from — oo to +oo along Im k > p and from +oo to —oo 
along Im k < —p. 

If ijj{k) is an entire function a{k), then (0) is the zero functional {a{k) 
need not be ultra analytic). 

ipo = j^dk a{k)ip{k) = (3) 

Equation (|^) tell us that ip{k) and ip(k) + a{k) represents one and the same 
functional ip{(p). 
If we take 

^Pik) = ^Sgilm k) (4) 
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then 

ipiv) = ^ J^dk Sg{Im k)ip{k) 

= - dk ip{k) - - / dk ip{k) 

^ J —CO Z: J +00 

dk <^{k) = l(v9) (5) 

The density given by (^) imphes that ip{if) is the unit functional. 

When we take ip{k) = for Im k > p (resp. for Im k < —p) the ultra 
distribution is the retarded (resp. advanced) functional. 

Feynman's causal function can also be represented as as ultra distribution. 
In fact, the propagator is: 

1 ^^cj>0 (6) 

kl — kP' — m? kQ — uj"^ 2uL!\kQ — UL! ko + uj' 

The causal Green function is such that the positive energy pole (at k^ = u) is 
propagated by a retarded potential and the negative energy pole propagates 
in an advanced way. To represent this Green function we take 

ipF{k) = — , for Im k > 

2uj k — uj 

tpFik) = —-r-^ — , for Imk <0 (7) 

If we interchange the poles, taking into account the senses of the path F, we 
get the anticausal Green function: 

1 1 

i^A{k) = -; , for Im k > 

2uj k + UJ 

1 1 

^Ak) = , for lmk<0 (8) 

2uj k — UJ 

Wheeler propagator (half causal and half anticausal) is then 

ilJw{k) = ^— ^ — - , for lmk>0 
2 k'^ — uj^ 

i^w{k) = -7;tT- — 2 ' lmk<0 
2 k'^ — 

, ,,, ISq(Imk) 
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We shall also need the Fourier transforms of the spaces ( and For a 
function (fi E (, the Fourier transform is: 

jTy, = ^(a;) = r^dk if{k)e-''''' (10) 

All functions if{x) are C°° and exp g|x|i^^<^(x) is bounded for any q and p. 
They form the space (. 

A natural definition for the Fourier transforms of the functionals ip E (' 

is: 

^^(^y.) = ^(^) = 27ii;{^) (11) 
For an ultra distribution represented by (|): 

2^T^P{(p) = iP{(p) = [ dk ij{k) dx ip{x)e''''' 
Jr J-oo 
+00 /• 

ikx 



i^i'f) = / dx (f{x) / dk ip{k)e 
J-oo Jr 

So that: 

/+00 
c/x ?/'(a;)(^(a;) (12) 
-00 

where 



^{x) = dk ^{k)e'^'' . (13) 

ip{k) is a density for (p{x) (g C)- Of course, if i^^k) = a{k) (entire analytic), 
ip{x) = 0. So that 'il}{x) is not altered if we add a{k) to ip{k). 

It should be clear that Eq. (|l^) is to be understood in the sense of dis- 
tributions. For example, if ip{k) is the density for the unit functional (Eq. 

(i): 

^(x) = - f dk Sg{Im k)e'^'' 

1 c+00 1 00 

= - / dk e''^" - - / dk e''^" 

Z J~oo L 



00 



^{x) = 2tt6{x) (14) 
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3 Pseudo-differential Operators 

Our aim is to work properly with some non-local pseudo-differential opera- 
tors. For this reason we are going to introduce the following operation on 
the functionals ip. 

Let us consider a function f{k) such that: f{k) is analytic in {k G 
C l\Im k\ > P} and f{k)/k^ is bounded continuous in {k G C/\Im k\ > /?}, 
j3 depends on /(/c), f3 e N. 

Then we define 

/ + 00 
-oo 

where: 

= J^dk f{k)ij{k)e''^ (16) 

and where the path F runs from — oo to +00 along Im k > (3 + p and from 
-|-oo to —00 along Im k < —j3 — p. 

We know that the functional ip{0) does not change when we add an 
arbitrary entire function a{k) to ip{k). However, such an addition in Eq. (|1^) 
gives rise to a new term 

A{x) = j^dk f{k)a{k)e'^'' (17) 

When f{z) is an entire function, A{x) = 0. For example for polynomial 
functions of the derivative operator. If f{z) = z~^, /(— '*^) is the inverse of 
the derivative (an integration). Equation ( p^ gives a primitive of ip^x). The 
additional term ( p!7| ) is: 

A{x) = [ dk ^e^fc^ = f dk ^e^^^ (18) 
Jr k Jl k 

where L is a loop around the origin. So that: 

A{x) = — 27rza(0) = arbitrary constant 
of course, an integration should give a primitive plus an arbitrary constant. 



/ 



dx 



■ip{x) 



0{x) 



(15) 
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Analogously, a double (iterated) integration, with f{z) = z ^, gives a 
primitive plus (p!?!): 



A{x) = / dk ^e'^^ = 7 + 



where 7 and 5 are arbitrary constants, as it should be. 

It is understandable that a more complex structure of singularities of f{z) 
gives rise to a more complicated A{x). 

To solve some pseudo-differential equations, we can work directly with the 
densities, the null functional having an arbitrary entire function as density, 
(an arbitrary constant has k~^a{k) as density, etc.) 

For a solution to a homogeneous equation: 

/(-^^)^(x) = (19) 

we write 

fik)^ik) = a{k) 

whose solution is: 

^(k) = r\k)a{k) (20) 
then replacing in Eq. (0) we obtain: 

i){x) = [ dk f~\k)a{k)e'^'' (21) 



In the case where the singularities of / ^ are concentrated on the real axis, 
the path V can be deformed to get: 



+00 



^{x)= dk f-\k + lO) - f-\k - lO) e"'^ (22) 
J —00 '- -' 

To solve the inhomogeneous equation: 

fi-^^M^) = Xix) (23) 
where x{x) is a given (generalized) function, we write 

fik)^Pik) = x{k) + a{k) 

m = r\k)x{k) + r\k)a{k) (24) 
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where x{x) = Jj-dk (x(^) + a{k))expikx = Jj-dk x{k)expikx. The last 
term in Eq. (plf) can be recognized as a general solution to the homogeneous 
equation (|T^). 

If in (p3D we choose x{x) = 5{x), we get the equation for the Green 
function: 

fi-'^Mx) = 6{x) (25) 
As a density for Dirac's function is \Sg{Im k) (cf. (|14D), we obtain from 

G{k) = ^f-\k)Sg{Im k) + r\k)a{k) (26) 

Eq. (|26D shows that, as usual, the Green function is defined up to a solution 
to the homogeneous equation. The different choices of a{k) determine the 
different prescriptions (advanced, retarded, Feynman, Wheeler, etc.) to avoid 
the singularities, when evaluating the Fourier transform. In the case where 
a = we have: 

G{k) = ]^r\k)Sg{Im k) + r\k)a{k) (27) 



It is interesting to find the Fourier transform of p7|): 

G{x) = 1 f dk r\k)Sg{Im k)e' 
2 Jr 



I. e. : 



Akx 



1 r+co , , 

G{x) = - dk (f-\k + iO) + f-\k - iO)) e' 

2 J — CO 

G{x) = g+ix)+g.ix) (28) 



where 



1 r+oo 

g± = - dkr\k±tO)e'''' (29) 

Z J— oo 

Let us now consider the composition of two operators /i and /2, to give 
fi-^i) = /i(-^i) ° For the densities we have /(fc) = h{k)f2{k) 

and practically nothing changes in the analysis we have made above. How- 
ever, if we apply the two operators in succession we have the following situ- 
ation: 
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For the equation 




(30) 



the square bracket is a solution of the homogeneous equation for /2. Accord- 
ingly: 



If we interchange /i and /2 in Eq. ( PPD , the solution (|HT|) changes is an obvious 
way. 

On the other hand, if /i and /2 are first composed to give / = /i o /a, 
the corresponding solution (|3l|) loses its last term. 

4 Some Examples 
4.1 Rational Functions 

We consider first a rational function (see also Ref. Ch. 2, §1.5), 




Now we have an inhomogeneous equation for /i, whose solution is 




(31) 




(32) 



n 



n 



P{k) = l[ik-k,) 



Q{k) = l[{k-k'^) 



(33) 



i=l 



{ki 7^ k'j for all i and j) 
The solution to the homogeneous equation: 




(34) 
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is given by ([T|) and 

■ip{x) - 

If we rewrite in the form: 



(35) 



(36) 



The solution 



is: 



■ip{x) 



dk 
dk 



Q{k) 
P{k) 
1 

m 



a(fc)e*'=^' + / dk 



P(A;) 



[g(A;)a(A;) + 6(A;)] e 



ikx 



(37) 



The square bracket in ( P?] ) is another arbitrary entire function. So that 

is in fact equivalent to (|35D. 
If now we write the equation: 



d , 



Q-^(-z^)^(a;) 



(38) 



The solution ( |HT| ) is: 



dk 



Q{k) 

p{ky 



a{k)e"''' + / dk Q{k)b{k)e 



ikx 



But the last integral is zero, as the product Qbexpikx is entire analytic. 

The three forms (P^), ( P3[ ) and ( |5B| ) are completely equivalent. The order 
of the operators P and Q is inmaterial. 

For the actual solution of (|35|), we take into account that 



1 



1 



E 



P(A;) U{k -h) Uk-k 



(39) 



where ai are appropriate constants (For the sake of simplicity we have as- 
sumed that all roots are simple). 
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The only singularities of the integrand in ( p5D are the poles at k = ki. 
The integral can then be evaluated by Cauchy's theorem: 



^(x) = [ dk ^^a{k)e'''' 



2'Ki ^ aiQ{ki)a{ki)e 

i=l 

n 



where j3i axe arbitrary constants. It is easy to check that 
(HD, (ll) and (H). 



(40) 

solves equations 



4.2 Power Functions 

For another (less simple) example, let us take the power function: 

where we introduce a cut along the negative real axis. 
For the homogeneous equation 

d 



(41) 



-%— ^(x) = 
ax 



we have: ip{k) = k "a{k) and 



ij{x) = I dk A;-"a(A;)e^'=^ 



(42) 



(43) 



The integrand is analytic, except on the (negative) real axis. We may write: 



?/'(x) 



dk 



{k + - {k - a{k)e 



ikx 



(44) 



But, (Ref. |23|, Ch. 1, §3.6) 



k-" + e-^'^^fcZ" 



{k + tOy 
[k-iO)-" = k-^ + e'^'^kZ" 



(45) 
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So that: 



+00 



iIj{x) = — 2isin(7ra;) / dk k_"a{k)e 
J —00 

The equation for a Green function is: 



ikx 



d \ 

-i— G{x) = 5{x) 
dx I 



I. e.: (cf. Eq. 



G{x) 



2 Jr 

1 r+00 



dk k °'Sg{Im k)e 



ikx 



dk 



(A; + zO)^° + (fc-iO)-" 



ikx 



For the Fourier transform we use Ref. |£3| (Ch. 2. §2.3) 



F(A;±iO)-" 



Via) 



Then: 



G{x) 



2T{a) 



(46) 



(47) 



(48) 



(49) 



To see the connection with the fractional derivative, we choose to add to 
), the solution (^) to the homogeneous equation with a{k) = —1/2. The 
new Green-function is the result of (j^T^ minus one half of 



G+{x) 



dk {k - iO)-"e 



a„ikx 



Via) 



-X 



a-l 



(50) 



Except of course for the exponential factor, Eq. ([50| ) coincides with the frac- 
tional derivative considered in Ref. [^(Ch. 1, §5.5). 

If instead of adding (^4]) with a = —1/2, we add the same function but 
with a = +1/2 we get: 



GJx) 



„a-l 



r(a) 



(51) 



which is also a fractional derivative operator. 



14 



5 Functions of the D'Alambertian 

In this section we are going to examine some non-local functions /(□) of 
□ = — 9^, where d' is the {d — l)-dimensional laplacian operator. 

For f{z) we will take an analytic function such as or {z + m^)", with 
a cut along the negative real axis, running from — cxd to zero or — m^, respec- 
tively. 

The ultra distributions (Eq. (H)) depend now on an [n — l)-dimensional 
vector k as parameter V'(^o) ~^ i^i^Q, k). 
The Fourier transform of Eq. (|TB|) gives: 

^{ko, k) = dko ^{ko, k)e''''^° (52) 

And, with the usual Fourier transform in the space of the parameters k, we 
obtain: 

^p{x) = iIj{xo,x = / d/c V(a;o, fc)e (53) 

J ~oo 

We define the operation /(□) on i(j{x) by: 

/(□)V^(x) =Jdk f{dl + Pmx,, k)e-"^' (54) 
And (Cf. Eq. (0)): 

f{dl + P)i^{x^, k) = dko f{-ko + Pmko, fc)e+^^°^° (55) 

We can now solve the homogeneous (n-dimensional) equation: 

f{a)^{x) = (56) 
According to (^4|), Eq. ( |56D implies: 

fidl + P)^Pixo,k) = Q 
And, due to (|55|) , we have for the densities: 

f{-kl + P)ij{ko,k)=a{kl + P) 
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k) = r\-kl + P)a(A;2 + (57) 

Where aik^ + k"^) is an arbitrary entire function of /cq, for any value of k. 
lYiom (p3) and (|57|) we get: 



'4j{xo, k)= dko r\-kl + l?)a{ko, A?)e*'=«^o (58) 



The analytic function of fco, / ^(— + presents a cut along the real axis, 
running from — oo to /cq = — t^, another twin cut from ko = +uj to +00 



(iu = +y k"^ + w?). 

The integration in (0) can then be taken along the real axis: 

^(xo, f ) = / dko A{f-')aiko, k)e'''^' (59) 
J —00 

where A is the discontinuity at the cuts: 

A(/-^) = f-\-{h + i^f + P) - r\~{ko - ^0f + P) (60) 

|a(/^"'^) = 0, for — uj < ko < uj , uj = +\[^ when m = o| 
For / = {-kl + + m^)" = {-kl + u;^)", we have: 

A(r^) = {-kl + u''-iQSg{ko))-''-{-kl + u^ + iQSg{ko))-'' 
= Sg{ko) [{-kl + uj^- «0)^" - {-kl +uj^ + 20)-"" 

And using (^: 

A(/-i) = 2i sin (7ra)5^A;o(-A;^ + cu^):" = 2z sin (7ra)5^A;o(A;2 - uj%^ (61) 



For a = 1, Eq. (|56D is Klein-Gordon equation. The weight function given 
by (pTf), has a pole for this value of a (Ref. |]23|, Ch. 3, §3.4). The residue is 
5{kl — u"^) = 6{k'^ — m^). For a ^ 1, Eq. (|6T|) gives A 6{k'^ — m^), which 
is the well-known invariant free wave solution. For a fractional a, the weight 
function gives a continuum with fc^ > m^. In this sense we can say that (|61|) 
represents free waves corresponding to a continuum of masses m . 

The free wave 6{k'^ — m^) concentrated on k"^ = rn? , changes for fractional a 
into (fc^ — m^)^°, which is spread from k"^ = to fc^ — > cxo (see also Ref. 
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For the fractional wave equation: 

= 

we can repeat the analysis made above. In this case: 



(62) 



A(/-^) 



i~kl + P- iQSgiko))-'' - i-kl + P + tOSg{ko)y 
Sgiko) U-kl - ^0)-" - i-kl + zO)-° 



A(/-i) = 2z sin {7ra)Sgko{kl)-'' (63) 

Again, for a — ^ 1, (|6^) is the usual wave equation and A(/~^) S{k'^), 
which is the elementary solution with support on the surface of the light- 
cone. When a is not an integer, Eq. (|63D implies that we can find free wave 
components for any positive mass squared. 

Let us now find a Green function for the fractional wave equation: 



□°G(x) 



dk {d^ + P)"G(xo, A?)e-'^"^ = Jdk (5(xo)e-'^"^ 
{d^, + PrG{xo,k) = 5{xo) 



And we have for the densities: 



Po + PTG{ko, k) = -Sglmko 



(64) 



Giko,k) 



k^Q + ky Sglmko 



G{xo,k) 



So that: 



G{xo,k) 



1 

2 J— oo 



dkn 



2 

' dko {-k^o + PySglmkoe''"'''^ 



[-k^ - + {-P + ioy 



(65) 
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The integrand in (^) does not depend on the sign of ko. This Wheeler Green 
function for the fractional wave equation (|6^) is an extension (for space-time) 
of the definition given in eq. (|^) . 

The two terms in the square bracket of (^5|) are, respectively, the causal 
and anticausal Green functions for the fractional wave equation. For a = 1 
they are the usual Feynman propagator and its complex conjugate. The 



Fourier transforms of those two terms can be found in Ref. Ch 3, §2.6. 
Using ( ^5]) we can write the square brackets as : 



[] = (A;2):- + e*™(A;2);"(A;2):" + 



so, 



1 r 
2 



-A;2-z0)-° + (-A;2 + i0)" 



(A;"):" + cos7ra(r) 



2\-a 



(66) 



Other Green functions can be found by adding to (|65|) solutions to the 
homogeneous equation. Adding together ([59| ) (with (0) and a = ±1/2) and 
(H) (with (HD), we get: 



dko (A;2):° + e 



^±^7r«Sgfeo^y^2^^-a 



ikQXo 



(67) 



The two terms of ( pSj) and the two Green functions of (|67|), coincide with the 
four Lorentz-invariant propagators for the fractional wave equation found in 
Ref. |11[. 

It is possible to repeat this procedure to find Green functions for the 
fractional Klein-Gordon equation: 



Instead of ( 

G{xo, k) 



j), we now have: 

1 /■+°° r 

dko (-A;^ + -iO)~" + (-P + + iO)" 



(68) 



2 J- 



Jkoxo 



(69) 

In this way we obtain the causal and anticausal propagators for the frac- 
tional Klein-Gordon equation. The Feynman function coincides with the 
propagator used in Ref. [0 to regularize the matrix elements of Quantum 
Electrodynamics. 
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6 Discussion 



The space ( of ultra analytic functions, defined in section 0, seems to be 
an appropriate base of test functions for the development of a framework 
adapted to the needs of those who want to work with non-local pseudo- 
differential operators. Its dual space C contains the propagators appearing 
in perturbative Quantum Field Theories. It is also flexible enough to allow 
for them a representation in terms of ultra distributions which is both, simple 
and general. 

The Fourier transformed spaces ( and translate the functions of C 
and (' into functions of coordinates where the pseudo-differential operators 
are supposed to act. In this way we can define non-local operations on 
(generalized) functions of the coordinates. These operators are extensions of 
the ordinary derivative operators, which are now particular cases in a wider 
scheme. 

We have written and solved homogeneous and inhomogeneous equations. 
In particular, we have found Green functions for some non-local operators. 
We have also shown the connection with the operators D"", defined and 



discussed in Ref . |]TT| , and with the fractional derivative of Ref . . In these 
references no attempt was made to solve the respective homogeneous equa- 
tions. We here exhibited the corresponding solutions. As a consequence, 
we are able to show that the free solutions of the wave equation, or the 
Klein-Gordon equation, with their sharp masses, are transformed when the 
equations are fractional, into a superposition of a continuum of masses with 
support in the interior of the light cone. Furthermore, the causal propaga- 
tors for the fractional equations, are seen to coincide with the analytically 
regularized propagators introduced in Ref. It is then possible to interpret 
this regularization method, as the matrix elements one would write for a non- 
local theory having no ultraviolet divergences. The usual infinities appear as 
poles for the local limit a ^ 1. 

We see then that with the chosen procedures we have convenient tools 
with which we can handle different pseudo-differential non-local problems. 
The scheme also provides a base for the quantization of non-local field theo- 
ries po 
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